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We derive a rigorous optical potential for electron-molecule scattering including the effects of 
nuclear dynamics by extending the common many-body Green's function approach to optical po- 
tentials beyond the fixed-nuclei limit for molecular targets. Our formalism treats the projectile 
electron and the nuclear motion of the target molecule on the same footing whereby the dynamical 
optical potential rigorously accounts for the complex many-body nature of the scattering target. 
One central result of the present work is that the common fixed-nuclei optical potential is a valid 
adiabatic approximation to the dynamical optical potential even when projectile and nuclear motion 
are (nonadiabatically) coupled as long as the scattering energy is well below the electronic excitation 
thresholds of the target. For extremely low projectile velocities, however, when the cross sections 
are most sensitive to the scattering potential, we expect the influences of the nuclear dynamics on 
the optical potential to become relevant. For these cases, a systematic way to improve the adiabatic 
approximation to the dynamical optical potential is presented that yields non-local operators with 
respect to the nuclear coordinates. 
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I. INTRODUCTION 

The concept of an optical potential as a physical entity that governs the scattering of a single particle by a composite 
target is an intuitively appealing phenomenological concept that goes back to the early days of nuclear physics. In 
principle, the scattering of a nonrclativistic, quantum-mechanical particle off an iV-particle target is a many-body 
problem and governed by the (N + l)-particle Schrodinger equation. In the so-called optical model [0, the elastic 
scattering problem is alternatively described by an effective single-particle Schrodinger (or Lippmann-Schwinger) 
equation. All effects of the interaction of the projectile particle with the target are contained in the so-called optical 
potential. In general, this optical potential has to be a very complicated object: It becomes a nonlocal operator 
because exchange and rearrangement of target particles have to be considered. An energy dependence has to 
account for possible excitations of the target and if inelastic scattering is energetically possible, the optical potential 
is nonhermitian in order to describe the loss of scattering amplitude into the inelastic channels. One major technical 
advantage of using optical potentials in numerical calculations is that the scattering problem can be separated from 
the many-body problem and the latter can be treated using bound-state techniques. 

The ability of many-body theory to derive exact, though not unique, optical potentials was a great success in the 
1950's |^-[||. The optical potential of Feshbach |l]-|3| follows from the most straightforward derivation by projection 
from the (N + l)-particle Schrodinger equation. An important alternative can be found in the self energy of the 
single-particle Green's function of traditional many-body theory . 

In contrast to the Green's function approach, the Feshbach optical potential is not easily amenable to perturbation 
theoretical approximations. Achieving a balanced and consistent treatment of ground- and excited-state correlation 
is thus difficult to achieve. Moreover, Mahaux and Sartor |?j have found in applications to nuclear matter that the 
nonlocality in the Feshbach optical potential is more complicated than in the Green's function optical potential. The 
numerical applications of Feshbach's theory on electron-molecule scattering reported so far, have been realised by 
projection of configuration interaction (CI) matrices |^,|9| but have rarely gone beyond the uncorrelated Hartree-Fock 
level for the description of the target's ground state. The inclusion of correlated, i. e., multi-configurational target 
wavefunctions in this context presents a delicate problem llOO]. 

On the other hand, the Green's function optical potential] ^["defined by the Dyson equation |l^] has turned out to 
be well-behaved and well-suited for numerical calculations. The Green's function approach to the optical potential has 
the advantage, that the self energy of the single-particle Green's function is easily accessible through the Feynman- 
Dyson perturbation expansion |12| . Thus a balanced treatment of ground-state and excited-state correlation is 
automatically achieved. Various approximation schemes for the self energies of molecules have been developed |lq] 
and applications to electron-molecule scattering have been reported in Refs. 16 - 20|. 
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All the previous calculations employing optical potentials have been done in the fixed-nuclei limit where the effects 
of nuclear motion are neglected altogether. A great variety of methods and numerical procedures can be used for 
solving the single-particle scattering equations in the fixed-nuclei limit like, e. g., the R- matrix pl|-|2"3|] , the Schwinger- 
variational [ p4| , the complex-Kohn f2aj2q] and the close-coupling [P7jp8l approaches. An often- used approximate 
treatment of the nuclear motion derives from the Born-Oppenheimer approximation. In this so-called adiabatic-nuclei 
approximation |29| , |30| , ^7| , fixed- nuclei S'-matrix elements for each configuration R of the atomic nuclei are averaged, 
weighted by the nuclear wavefunctions of the isolated target molecule. This approximation, however, breaks down 
for low projectile energies and, in particular, for resonant vibrational excitation (see e. g. Finally, in situations 

where the electronic projectiles move with velocities that are comparable to the typical velocities of nuclear motion, 
it becomes indispensable to treat the coupled motion correctly. 

The most straightforward way of accounting for the full, nonadiabatic coupling of the nuclear motion to the 
motion of the projectile electron is realised in the close-coupling formalism p^ , |2^ ]. In this formalism, the scattering 
wavefunction of the full problem (i. e. the degrees of freedom of N + 1 electrons and the nuclear degrees of freedom 
appear as variables) is expanded in a finite set of basis functions that are chosen according to appropriate exact 
or approximate constants of motion. Usually it is not feasible to include excited electronic states of the target in 
the close-coupling expansion. These states, however, are necessary to describe polarisation effects correctly. In the 
literature, so far model potentials void of rigorous justification have been used to account for this defect |B7LB3 -Bi). 
The close-coupling formalism further has the problem that the numerical calculations are difficult to converge for very 



low projectile energies |35 



An alternative to the one-step close-coupling expansion for treating the nuclear dynamics in electron-molecule 
scattering is a two-step procedure employing the so-called projection-operator (PO) formalism ]3^-|3^|. Based on 
the assumption of a resonance in fixed-nuclei electron scattering, the PO approach can be used to treat the nuclear 
dynamics going beyond the adiabatic-nuclei picture. In the PO formalism the electronic scattering problem is separated 
from the nuclear dynamics by Feshbach projection such that only some resonance parameters, which can be taken 
from fixed-nuclei scattering calculations, determine the equations for the nuclear dynamics that have to be solved. 
The coupling to the continuum of electronic scattering wavefunctions enters the nuclear dynamics via a complex, 
nonlocal potential. Different levels of approximation for the nonadiabatic coupling of electronic and nuclear motion 
are reached by local and nonlocal models for the complex potential The PO method has been very successful 
in explaining many features of resonant electron-molecule scattering and it gives a very intuitive picture for the 
mechanism of many physical processes like, e. g., resonant vibrational excitation j37]]. On the other hand it is per se 
not an exact approach because only resonant contributions to the cross sections can be treated. There are always 
non-resonant contributions (also called background contributions) that cannot be evaluated in the PO formalism but 
may actually dominate the spectra in off-resonance scattering. Sometimes no dominating resonance is present at 
all, and sometimes numerous resonances contribute which are difficult to identify individually. For these cases the 
separation of nuclear and electronic motion ultimately fails. Also the nonlocality of the complex potential usually 
becomes very complicated and difficult to model for very slow scattering electrons: Details of the electron-molecule 
interaction and the dynamical coupling of electronic, vibrational, and rotational degrees of freedom become important 
when the time scales of projectile and nuclear motion in the target are comparable. As an example we may mention 
the scattering by a polar molecule with a supercritical or near-to-critical dipole moment. As it is well-known, a 
supercritical, nonrotating dipole (D > 1.625 Debye) can bind an excess electron but this bound state is weakened or 
even disappears if the dipole is allowed to rotate 39 1. It is, however, very difficult to include the correct coupling 
of rotational degrees of freedom in the nonlocal complex potential model jlO| and thus a full dynamical calculation 
including all relevant degrees of freedom may become necessary in order to describe the physics of very-low energy 
scattering by polar molecules correctly. 

To our knowledge, no one has so far given a rigorous derivation of an optical potential for the scattering electron 
that takes full account of the nuclear motion. To close this gap is the main objective of the present communication. In 
this paper we will present an optical potential that is based on the usual Green's function optical potential but allows 
to treat projectile and nuclear motion on the same footing. This so-called dynamical optical potential is given by 
the self energy of the dynamical Green's function, an extension of the usual single-particle Green's function including 
the nuclear kinetic- energy operator. The dynamical optical potential is an optical potential for electron-molecule 
scattering that is clastic with respect to electronic excitations of the target but inelastic with respect to the vibrations 
or rotations of the target. The only assumption that is used throughout this paper is that the Born-Oppenheimer 
approximation is valid for the target's ground state. Otherwise the dynamical optical potential is exact. It takes 
full account of the many-body effects of electron scattering including exchange and polarisation as well as the full 
nonadiabatic coupling of the projectile motion to the nuclear motion of the target. The expressions we derive for the 
exact optical potential may provide guidance for the choice of model potentials on the one hand, and on the other 
hand are amenable to ab-initio approximations and can be calculated with standard quantum-chemical methods. 

This paper is organised as follows: In the next chapter we will start with reviewing the definition of the traditional 
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single-particle Green's function of many-body theory. Applied only to the electronic degrees of freedom of a mole cule 
at fixed nuclear configuration we also call this propagator the electronic Green's function. In the subsequent Sec. [I B 
we then define the dynamical Green's function, which allows for the treatment of the full molecular dynamics including 
the nuclear degrees of freedom as dynamical variables. We also define the closely related (nuclear) inelastic Green's 
function that shows the close relation to the recently developed theory of electronically inelastic Green's functions. 
Sec. |j] is closed by showing how the electronic propagator can be obtained in a suitable limit from the dynamical 
Green's function. In the third chapter an algebraic derivation of Dyson's equation for the electronic Green's function 
is reviewed before the dynamic Green's function is shown to fulfil a Dyson equation and the dynamical self energy 
is defined. Chapter IV treats the relation to electron-molecule scattering. The dynamical self energy is shown to be 
closely related to the optical potential for the coupled motion of the projectile electron and the atomic nuclei of the 
target molecule in electronically elastic scattering. Expressions for the scattering S and T matrices are shown for 
inelastic processes with respect to the nuclear coordinates including vibrational or rotational excitation, associative 
detachment, and dissociative attachment. In Sec. IV B[ an effective Schrodinger equation using the dynamical optical 
potential is d iscus sed. A direct derivation for the effective Schrodinger equation based on projection operators is 
given in Sec. 



rvc 



which yields additional insight into the relation of our Green's function approach to Feshbach's 
theory. The well-known stat ic-exch ang e and th e fixed-nuclei sc atterin g equations are identified as approximations to 
our exact equations in Sees. [VP 1 and IV D 2 . The rest of Sec. IV D is devoted to the discussion of different possible 
approximations to the dynamical optical potential, in particular concerning the treatment of the nuclear dynamics 
in virtual electronic excitations of the target. We also discuss how to compute the dynamical optical potential using 
standard ab initio quantum chemistry methods. 



II. THE DYNAMICAL SINGLE-PARTICLE GREEN'S FUNCTION 

Traditional many-body theory [ p"2[ deals with systems of indistinguishable particles. Applications to molecules 
usually treat only the electronic degrees of freedom in the many-body formalism and freeze the nuclear configuration. 
After briefly reviewing the definition of this purely electronic single-particle Green's function we will define the 
dynamical Green's function, which allows us to treat the nuclear degrees of freedom as dynamical variables. 



A. The electronic Green's function 

The usual single-particle Green's function G{rt, r't') of non-relativistic many-body theory is defined by |Q 

iG(rt, r't') = «| ^(r, t)^(r', t') |<) 0(t - t') 

-{^\^{r',t')^{r,t)\^)9{t'-t). (1) 

Here, 6(t) is Heavyside's step function, l^^} denotes the exact (correlated) ground state of the TV-particle system, 

and ift (r,t) is the Heisenberg field operator, which creates a particle at position r and time t. Although we consider 
fermionic particles, electrons in particular, we suppress the notation of spin indices for simplicity (r may be considered 
as spin-space variable) . In applications of many-body theory to molecular physics one usually assumes an adiabatic 
decoupling of electronic and nuclear motion and thus only treats the electronic degrees of freedom explicitly. G(rt, r't') 
is then the Green's function in a system of electrons moving in the static potential of the nuclei, fixed at the coordinates 
R. In the following, we will call G(rt, r't') the (purely) electronic Green's function in contrast to the later defined 
dynamical Green's function, which will allow to treat non-adiabatic coupling of electronic and nuclear motion. 
The Green's function decomposes into two parts, the so-called particle propagator G + and hole propagator G~: 



G(rt,r't') =G + (rt,r't') + G-(rt,r't'). (2) 

Noting that the (purely electronic) Heisenberg operator if) (r,t) evolves in time with the electronic Hamiltonian H c 

^\r ) t)=e iH ^(r)e- iH ^ ) (3) 

we can rewrite the particle and the hole propagator in the following way: 

iG+(rt,r't') = {^^\i;(r)e- l[H ^ Va ^ t - t '^\r')\^)0(t-t'), (4) 
iG-(rt,r't') = -^^^(ry-^^-^^'mr) \^%)9{t' -t). (5) 
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Here, we made use of the fact that the electronic ground state \^q) is an eigenstate of the electronic Hamiltonian 
H e with the energy of the electronic ground state Vo(R): 

H C \^) = V (R)K). (6) 

We explicitly indicate here the dependence of Vq on the nuclear coordinates R for later use. Of course, H c and 
depend parametrically on R as well. 

For convenience, we will later discuss the Green's functions mainly in the frequency instead of the time domain and 
using an orbital instead of a position space representation. The corresponding expression for the electronic Green's 
function G pg (ui), which is equivalent to the definition (|l|), reads 

Here, the Fourier transformation into the frequency domain is defined by F(u>) = f^^dft — t') e IW ( t_ * ^F(t — t'). 
The positive infinitesimal 77 remains from a convergence factor that has to be introduced to assure the convergence 
of the Fourier transformation. The sign of the term irj determines the time ordering by theta functions in Eqs. (Q) 
and (|5|) jl^] . The transformation into the orbital representation is defined by the transformation of the creation and 
destruction operators of second quantisation. For convenience we assume that {<f p (i')} p is a discrete but complete set 
of square-integrable and normalised basis functions for the single- particle Hilbert space. The creation operator for an 
electron in the orbital <p P (r) is defined by 
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dr<p p {r)^(r). (8) 



Like usual, the orthonormality of the orbitals implies the canonical anti-commutation relations {aj, , a q } — S pq , etc. 

The orbital set has been chosen discrete just for conceptual clearness and simple notation. The generalisation to 
continuous sets does not pose principal problems. Convenient realisations of orbitals may be found, e. g., in the 
Hartree-Fock orbitals of the target (together with a convenient discretisation of the continuum) or in the momentum 
space representation (discretised by placing the system in a finite box and employing periodic boundary conditions). 



B. The dynamical Green's function 

If we now want to couple the electronic motion dynamically to the nuclear motion we have to replace the electronic 
Hamiltonian H c by the full molecular Hamiltonian 

H = T n + H c . (9) 

It is given by the sum of the nuclear kinetic energy operator T n and the electronic Hamiltonian H e . By convention, 
H e contains the internuclear repulsion in addition to the electronic kinetic energy as well as the electron-nucleus 
Coulomb attraction and the electron-electron repulsion. The nuclear kinetic energy operator T n generally includes 
all vibrational, rotational and translational degrees of freedom for the target molecule. For convenience we assume 
that the translational motion has been separated off. In this case rotation, vibration, and dissociation are the 
remaining degrees of freedom for the nuclear motion. In any case, T n is a differential operator that acts on the nuclear 
coordinate R. Thus T n does not commute with H c or Vq(R). However, we want to assume in the following that 
the Born-Oppcnhcimcr approximation is valid for the electronic ground state of the molecular target, which implies 
that the commutator of T u with the electronic ground state [T n , \^q)] is negligible. Physically, this means that the 
electronic configuration in \^q) depends only weakly on the nuclear coordinates. 
Let us introduce the nuclear Hamiltonian 

H n = T n + V Q (R) (10) 
and let 4> k (R) be an eigenfunction with eigenvalue Eok- 

H n ct> k {R) = E ok 4> k {R). (11) 
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In the following we will also use the abstract Dirac notation \<pk) from time to time for the nuclear eigenstates. The 
wavefunction 4>k{R) then emerges as the coordinate representation of the state \4>k) by (f>k{R) = {R\4>k)- 

Under the assumption that T n commutes with it is easily seen that the product wavefunction \^Q)ipk(R) 

becomes a molecular eigenstate with energy Eq^: 

H = E ok \^)MR)- (12) 

This Born-Oppenheimer picture is usually an adequate description for closed-shell molecules in the electronic ground 
state. The negligibility of the commutato r \T n , \^q)} is the only fundamental assumption of the theory we will 



develop in the following. At the end of Sec. IV C we will return to discuss the validity of this approximation and show 
how the theory can be generalised. At this point we want to stress that we only assume the Born-Oppenheimer picture 
to be valid for the nuclear dynamics in the electronic ground state of the molecule. We neither assume adiabaticity 
for the scattering electron nor for the electronically excited target states! 

We now define the dynamical Green's function Q(rt,r't') as the sum of the two components for particle and 
hole propagation: 

g(rt,r't') = g+(rt,r't') +g-{rt,r't'). (13) 

The particle part Q + is defined by 

ig + {rt, r't') = (tf^l ?A(r)e- ,;(ff - Boo)(t -*' VV) 1*^} 0(t - t'). (14) 

This propagator resembles the standard particle propagator (Q) where the electronic Hamiltonian H c appears instead 
of the total Hamiltonian H and the electronic ground state energy Vq(R) replaces the molecular ground state energy 
Eqq. Note that H now includes the nuclear kinetic energy T n , which is a diffe rential operator on the nuclear coordinates 



R. Thus, also g is an operator on R. As will be shown explicitly in Sec. IV D 4, the dynamical Green's function g 



may also be understood as a nonlocal integral operator in the nuclear coordinates: 

(R\ g{rt 7 r't') \<f>) = J dR'g(rRt, r'R't') <f>(R'). (15) 

For the sake of simple notation, however, we will not indicate the Independence of g explicitly. 

The hole part g~ is defined with a modified Hamiltonian because this will be helpful later for the formulation of 
the Dyson equation: 

ig-(rt, r't') = ^\^\r'Y^- Eaa - 2 ^~ E ^ t - t '^{r) 6{t' - t). (16) 

As long as we describe scattering or attachment processes, only g + has physical significance and the propagation 
with the "wrong" Hamiltonian H — Eqq — 2(H n — Eqo) does not enter the physics. If, however, we want to study 
ionisation where the hole part carries physical significance, we have to modify the definition of the dynamical Green's 
function g and allow g~ to propagate with the correct Hamiltonian H — Eqq. In this case we may modify g + in order 
to obtain a well-behaved Dyson equation. In the following we will keep to the choice ([u]) and ( |l6| ) because we are 
mainly interested in the scattering problem in this communication. 

Note that the particle part g + can still be expressed by an expectation value of a product of Heisenberg field 
operators: 

ig+(rt,r't') = e - l (^- £ oo)^JV| ^ rt )^(rt') \^)e l{H ~- Eaa)t ' 6{t ~ t'). (17) 

In contrast to the definition of the purely electronic single-particle Green's function ([!]), the field operators ip(rt) and 
tp' (rt') now evolve in time with the full molecular Hamiltonian H = T n + H c and thus feel the effect of nuclear motion 
in the target. The hole part g~ cannot be expressed in such a simple way as g + . 

We conclude this subsection with the expression the dynamical Green's function in the frequency domain and orbital 
representation: 

6 M M = «|a ^_ g+ 1 EM+ . |) 4|<) 

This expression is equivalent to the definition of the dynamical Green's function in the time domain, ( |l3| ) to (p^). 
The transformation to frequency and orbital representation is defined as before [see Eq. (0)] . 
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C. The inelastic Green's function 



The dynamical Green's function Q presents, like the purely electronic propagator G, a matrix (or kernel of an 
integral operator) in time and in the coordinates of a single electron. Moreover, Q is an operator in the nuclear 
coordinates R. For later use and for comparison with the (electronically) inelastic theory of reference , we define 
the inelastic Green's function C?l m > n ] as the matrix element of Q with respect to the nuclear wavefunctions <p n (R): 

glmri^tyt') = (ct> m \G{rt,r't')\cl> n ). (19) 

This inelastic Green's function relates directly to the S-matrix of rotationally or vibrationally inelastic scattering like 
we will see later. 

The relation to the inelastic Green's function studied by Cederbaum becomes apparent when we evaluate the 
particle part of g'™ 1 '™'. We introduce the symbolic notation \0k)) for the molecular state \^Q)4>k{R) of Eq. (|T^) , 
which describes a (ro-) vibrational excitation with quantum number k in the electronic ground state. From Eqs. Qll| ) 
and (Il7|) it follows that we can write the particle part of the inelastic Green's function as 



i0K»]+( rt ,r'f') = «0m| ip{rt)^(rt') \o n )) e ^K E o m -Eo )t^(E 0n -E 00 )t'] d ^ _ ^ ^ 

where the double bracket notation ((-|-)) indicates integration over both electronic and nuclear degrees of freedom. 
Expression ( |20| ) is completely analogous to the definition of Rcf. 41 1. While in this Reference inelastic scattering 
between different electronic excitations of the target is considered and the nuclear degrees of freedom are neglected, 
the present study focuses on inelastic processes with respect to the vibrations and rotations in the same adiabatic 
electronic configuration. In contrast to Ref. jyj where the projectile particle, for instance a positron, is distinguishable 
from the target electrons, our projectile is an electron and we fully account for the indistinguishable nature of the 
projectile and target electrons in the present approach. For this reason we need the hole part Q~ , which, admittedly, 
takes on a less intuitive form than the particle part Q + but will prove very useful later for the formulation of the 
Dyson equation. 

We now want to show that the dynamical Green's function Q boils down to the usual electronic Green's function G 
when the nuclear dynamics becomes unimportant. Formally the transition from Q to G is achieved by assuming that 
the Hamiltonian of nuclear motion H n = T n + Vo(R) can be replaced by its lowest eigenvalue Eqo. When we express 
the nuclear kinetic energy T n in the molecular Hamiltonian H — T n + H c by H n — Vq (R) and replace H n by Eqq , the 
parts of the dynamical Green's function Q + and Q~ in Eqs. (|lj) and ( |l6[) reduce to the conventional expressions (^) 
and (^|) for the electronic Green's function G + and G~ . Also in the inelastic formalism the transition to the purely 
electronic Green's function can be achieved. The elastic channel component Q\°' ^ can be expressed through the purely 
electronic Green's function under the assumption that the nuclear kinetic energy T n commutes with the electronic 
Hamiltonian H c : 

0[°'°](rt,rV) lT »'JH^° J dRcp*(R)G(rt,r't')MR)- (21) 
For the inelastic components we find an expression that reminds us of the adiabatic-nuclei approximation: 

g[n,m] I T » j dR e - iiEom - Eoo)t (P* n (R) G(rt, r't') cj> m (R)e^ Eo '>- Eoo)t . (22) 

The phase factors account for the different reference energies of the channel states. 



III. DYSON EQUATION AND DYNAMICAL SELF ENERGY 



In this section we will derive Dyson's equation for the dynamical Green's function and we will show that the 
dynamical self energy can easily be expressed as a generalisation of the purely electronic self energy. Before turning 
to the dynamical Green's function, however, we will briefly review the algebraic derivation of Dyson's equation for the 
traditional single-particle Green's function. The Dyson equation for the dynamical Green's function will then follow 
easily as a generalisation of this procedure. 



G 



A. The Dyson equation for the electronic Green's function 



The following derivation of the Dyson equation is close in spirit to Refs. [|42| and 43 1. In the present form it has 
already been given in Ref. |44| . For convenience we use the representation of the Green's functions in the frequency 
domain and orbital representation in this Chapter as defined by Eqs. (0) and ([l8|). 

The definition (]?]) of the electronic Green's function may be written in a more compact form, combining the particle 
and the hole propagator. To this end we introduce the following composite vectors 

and the matrix 

H -{ Vo(R)-H s )- ^ 

We can rewrite Eq. (Q) by formally expressing the single-particle Green's function as a matrix element of the 
resolvent of H with respect to the states 

G m {u) = (Y p \—^\Y q ). (25) 
Here, we omitted the ir\ terms to avoid clumsy notation. The time-ordering controlled by these t erms i s usually not 



important but may become relevant in time-dependent formulations of scattering theory (see Sec. IV A). 

The scalar product of the composite vectors \ Y p ) now includes the matrix- vector product where bras always move to 
the left and kets move to the right. With respect to this canonical scalar product, the vectors \Y p ) fulfil the following 
orthonormality relation: 

(Y p \Y q ) = «| 0pffl t |^) + (^| a t ap |^) = §pq . 

When the iV-electron wavefunction \^q) is approximated by a Slater determinant, the vector \Y p ) contains either 
a single-hole state, if p refers to an unoccupied orbital, or a single-particle state, if not. In the more general case of 
a correlated wavefunction jv^o ), the vectors \Y p ) mix N — 1 and N + 1-electron states. Note that in either case and 
for arbitrary single-particle indices p, the vectors \Y p ) may be seen as orthogonal states that span a linear space with 
the same dimension as (or isomorphous to) the single-particle Hilbert space. We will call the space spanned by the 
vectors {1^}};) the primary space. Since the primary states \Y P ) and their linear combinations are not, in general, 
exact eigenstates of the electronic Hamiltonian H c , they couple to the higher excitations like two-particle-one- hole 
excitations, one-particle-two-hole excitations, etc. The space of higher excitations will be called the secondary space 
fj5"f . The concept of single-hole, two-particle-one-hole excitations, etc. is, of course, only adequate when the target 
state is dominated by a single configuration. In general, \^q) will be represented by a correlated wavefunction 

and thus the primary and the secondary space cannot be easily expressed by single-configuration basis states. For 
an explicit construction of a basis for the secondary space that allows for convenient approximations to the Green's 
function or the self energy in the framework of the so-called intermediate state representations see Refs. ^,|^|. In 
the following we assume that we have a basis of the composite space { \Qj)}j that consists of the primary states 
{|Yp)} p augmented by any suitable basis for the secondary space. The particular choice of the basis for the secondary 
space does not matter. 

We may now proceed to derive the Dyson equation and define the self energy setting out from Eq. (^5|). The basis 
{ \Qj)}j of the composite space defines a basis set representation H_ of the matrix Hamiltonian H: 

mu = (Qi\H\Qj). (26) 
By virtue of the subdivision of the basis set into two parts, the matrix f£ is structured into blocks: 

The block index a refers to primary states VY P )} P and b to the rest of the basis { \Qj)}j- The upper left block of 
this matrix is readily evaluated using Eqs. (|23) and (|3): 
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H 



(Y p \H\Y q ) = (^\{a p , [J?, a t]}K>. (28) 



It is easily seen, now, that the electronic Green's function of Eq. (g5|) can be understood as the upper left block of 
an inverse matrix: 

\ uil — H I 

\ = aa 

The proof follows by using the completeness of the basis { \Qj)}j and the orthonormality of the primary states 
{|^p)}p- By simple matrix partitioning, the inverse matrix of the Green's function G(u>), can be expressed by 

CM -1 = wl -H -H , -^—H (30) 

In order to define the self energy and to make contact to the usual form of Dyson's equation, we introduce a 
perturbation theoretic separation of the electronic Hamiltonian H e into a zeroth order, an interaction part, and the 
internuclear repulsion term Un(R), which does not act as an operator on the electronic coordinates: 

H e = H efl +H e>1 + U N (R). (31) 

The zeroth order part 7J e ,o is a one-particle operator of the form 

//,,, •.-"."/• ( 32 ) 

A convenient choice may be realized by the electronic kinetic energy or the Hartree-Fock operator. We do not demand 
the operator H c ^ to be diagonal in the given orbital basis but a diagonalising choice of the single-particle basis is 
always possible. E. g. if -ff c .o is chosen to represent the kinetic energy, the matrix e with the matrix elements [el . . = 
becomes diagonal in the momentum representation. In atomic units, the diagonal element is then simply given by 
e p = p 2 /2. In practical applications, H e $ is often chosen the Hartree-Fock operator. 

On the assumption that the electronic ground-state wavefunction \^q) and energy Vq(R) both possess well-behaved 
perturbation expansions, the zeroth order of the matrix H_ can be easily evaluated from Eq. ( ^jj ) to give the matrix 
of the zeroth-order orbital energies: 

£ (0) = e. (33) 
It is also easy to see that the zeroth orders of the off-diagonal blocks H , and H, vanish and hence the zeroth-order 

J ° =ab =ba 



electronic Green's function G' ' is found from Eq. (pG) to yield 



G<°>H = -^— . (34) 

— COL — £ 



The electronic self energy is defined by 



1 



Y.(iu)=H -e + H , ^H, (35) 



and Eq. (pG) remains 



= u i-g-m*y (36) 

This is equivalent to the common form of Dyson's equation: 

G(u) = G (0) (w) + G (0) (w)S(w)G(w) . (37) 
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B. The Dyson equation for the dynamical Green's function 



We will show in the following that a Dyson equation for the dynamical Green's function Q can be derived in an 
analogous manner, using the concepts of the previous section. We will discuss in particular the dependence on the 
nuclear coordinates, which becomes important because of the introduction of the nuclear kinetic energy operator. 



Just like the electronic Green's function, the dynamical Green's function Q of Eq. (18) can be written in a compact 
form combining the ionisation and attachment parts. Using the composite states |Y^,) of ( |23| ) and the electronic matrix 
Hamiltonian H of (p4|), we can rewrite the expression (|l8|) for the dynamical Green's function to give 

0p» = ( y p\ TT~^ ^H^>- ( 38 ) 

We again dropped the ir\ terms here for brevity. Also the 2x2 unit matrix that allows u>, H n , and E aQ to be applied 
to the two-component vector \Y P ) is omitted. The particularly simple and compact form of this expression can be 
seen as a preliminary justification for the particular choice of Q~ in the definition (|l6|). Dyson's equation can now 
easily be derived in an analogous way to the last section. 

Comparing with the expression (|25| ) for the electronic Green's function, we see that only the term — H n + Eqq is 
additionally present in the dynamical Green's function of Eq. (|38|). The molecular ground state energy Eqq is just 
a constant and defines the zero point of the oj scale. The "nuclear" Hamiltonian H n = T n + Vq(R) introduces the 
nuclear kinetic energy operator T n , which is a differential operator on the nuclear coordinates R. It is therefore useful 
to consider the H-dependence of the quantities used in the derivation of the Dyson equation. 

The matrix operator H contains the electronic Hamiltonian H c and the ground state energy Vq(R). Both quantities 
depend on the nuclear coordinates R and thus care has to be taken because H and H n do not commute. For the 
moment we want to assume that the electronic orbital basis that defines the electronic creation operators al is 
independent of the nuclear coordinates R. This is the case, e. g. for the momentum or position space representation. 
This restri ction is just convenient for the derivation of the Dyson equation but not essential and we will see later 



(Sec. |IV C| ) that it can be lifted as long as the chosen single-particle basis for the electrons is complete. The R 
dependence of the composite s tate s |T^,) of Eq. (p3) now derives entirely from the electronic ground state \^q). Like 



we already mentioned in Sec. II B, we assume that the Born-Oppenheimer approximation is valid for the electronic 
ground state and thus the nuclear kinetic energy T n commutes with | ) . Consequently we may assume that the 
nuclear Hamiltonian H n commutes with the composite states |Yp). We also assume that the basis states \Qj) of the 
secondary space are conveniently chosen, such that the commutator with T n or H n can be neglected. An appropriate 



choice is always possible 47 . The nuclear Hamiltonian H n can consequently be pulled out of any matrix clement 
involving the basis states of the composite electronic space and thus the matrix representation of H n in this basis is 
proportional to the unit matrix: 

(Qi\H n \Q J )=S I jH n . (39) 

In analogy to Eq. ( p9| ) we can understand the dynamical Green's function Q as the upper left block of an inverse 
matrix: 

£M=h rr \ s i • ( 4 °) 



(u-H n + E 00 )l-H 



aa 



The elements of this matrix are now operators acting on the nuclear coordinates. Matrix partitioning can be applied 
and is completely analogous to the case of the electronic Green's function. 

Owing to the partitioning of the electronic Hamiltonian H c into a zeroth order and interaction part, we find 
for the zeroth order dynamical Green's function 

£ (0) M = 7 v , (41) 

where e is again the matrix of zeroth-order single-electron energies. In the simple case where the electronic kinetic 
energy is chosen as the zeroth-order electronic Hamiltonian -ff c ,0; the- matrix e is diago nal in the momentum repre- 



sentation and independent of the nuclear coordinates R as mentioned in the Sec. Ill A. The zeroth order dynamical 
Green's function now is also diagonal and describes the simultaneous motion of a free electron and the vibrations 
(or rotations) of the isolated target molecule. 

In the more general case where the matrix e depends on R (if, e. g. -£/c,o is chosen the independent Hartree-Fock 
operator), the zeroth-order propagator may be understood to describe the coupled motion of a vibrating (or 







rotating) molecule in the electronic ground state and a single electron that moves under the influence of il-dependent 
mean fields. Since the zeroth-order electronic Hamiltonian H c _q is a one-particle operator, the motion of the scattering 
electron is decoupled from the target electrons. The effective Hamiltonian for the coupled projectile and nuclear motion 
is thus independent of the scattering energy and simply given by H n + iJ c .o — Eqq. The zeroth-order Green's function 
Qt ' is the resolvent of this effective Hamiltonian. 
For the full dynamical Green's function we find 

^ M = {.-H n + E^)l = -e = -A^y (42) 
where the dynamical self energy A is defined by 

A(ui) = H -e + H . =— JT . (43) 

=aa = =-"( W -i/ n + S o)l-^= fcQ 

In comparison to the zeroth order case, the dynamical self energy A(u>) now accounts for the many-body nature of the 
molecular target. We will demonstrate in the next section that A(u>) indeed is an optical potential for (electronically 
clastic) electron-molecule scattering. It is remarkable that the nuclear Hamiltonian H n appears in the self energy 
A(uj) together with the uj dependence and thus introduces derivatives with respect to the nuclear coordinates. Like 
we shall see in detail later, the dynamical self energy thus becomes a non-local operator not only in the electronic but 



also in the nuclear coordinates. In Sec. IV D we will discuss the meaning of the nuclear Hamiltonian appearing in the 



dynamical self energy and possible approximations. 

In a very formal manner, A(uj) can be expressed by the usual electronic self energy S(w): 

A(u)=UuJ-H n + E 00 ). (44) 

Before we turn to the discussion of the scattering S'-matrix and the optical potential in the next section, we would 
like to present the Dyson equation in terms of the inelastic Green's function £/[ m ' n l of Eq. (|l9|). Carrying on the 
analogy to the inelastic formalism of Ref. |4l| ] will help us in the next section to identify the dynamical self energy 
with the optical potential. The inelastic formalism also yields a slightly different view of the Dyson equation as 
the introduction of nuclear cigcnfunctions <j) n (R) formally leads to an equal treatment of the electronic and nuclear 
degrees of freedom, at the cost of having to introduce another sort of matrix indices. 

Instead of using the purely electronic basis { \Qj)}j, we introduce the product basis { \Qj) <S> \4> n )}j,n, where \<j) n ) 
is the abstract notation for the nuclear wavefunction cf> n (R) . The expression (|l9|) for the inelastic Green's function in 
the frequency domain and orbital representation then reads 

Q\™>"\ U ) = (d> m \ (Y p \ g-J- —\Y q ) \cj> n ). (45) 

The bracket of the nuclear states (4> m \ ■ ■ ■ \<fin) implies integration over the nuclear coordinates. The derivation of the 
Dyson equation can now be redone using the matrix representation defined by the product basis. In the following we 
denote matrices in the labels of the nuclear states n, m by boldface letters and matrices in the electronic state labels 
/, J, p, q with double underbars. The zeroth order inelastic Green's function becomes 

g' ' ( " ) = (w + ^)l- g 4-I - (46) 

The matrix of the nuclear Hamiltonian H n is diagonal due to the particular choice of the nuclear wavefunction basis: 

[Hj m ^ =S mn E 0n . (47) 

The full inelastic Green's function reads 

^ ] = (c + ^I-inl-l-^)' (48) 

where the inelastic self energy A(lo) is defined by 

A(uj)=H -e + H u ^H, . (49) 

= =aa = =afc (^ + ^ 00 )l : - J ff n l-^ b& =ba 
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Of course, matrix inversion and multiplication now have to include both sorts of matrix indices. The Dyson equation 
for the inelastic Green's function can be written as 

g(u) = g(°) ( u ) + ( u ) £fw) g{oj) . (50) 

Note that the inelastic self energy A is merely another representation of the dynamical self energy A. The relation is 
given by 

A™< n \u) = {<t> m \A m {u>)\<t> n ). (51) 



IV. APPLICATION TO ELECTRON-MOLECULE SCATTERING 

In this chapter we discuss the relation of the dynamical Green's function to the process of electron-molecule 
scattering. We introduce the S- and T-matrices of electronically elastic but vibrationally or rotationally inelastic 
scattering and show how they can be calculated using the dynamical Green's function formalism developed in the last 
sections. We present an effective Schrodinger equation and discuss the role of the dynamical self energy as an optical 
potential for the scattering process. An direct derivation of the effective Schrodinger equation is given and possible 
approximations to the exact equations are discussed. 

A. The S- and T-matrices of inelastic scattering 

The relation of the dynamical Green's function to the S- and T-matrices of scattering theory is very similar to 
elastic electronic scattering off atoms or rigid molecules SB] and completely analogous to the electronically inelastic 
case jy} . We therefore sketch the derivation only very briefly. 

The S'-matrix for the processes we consider reads 

S{p'm^pn) = {{V;r\% + )), (52) 

where the stationary scattering states I 1 ?™ 1 * 1 )) with incoming (— ) or outgoing (+) boundary conditions are defined 
by 

I*"*)) = lmi e -*(«p+*>»-tf)* t |**) |^), (53) 

y t— >=foo 

and the double bracket notation again implies integration over electronic and nuclear degrees of freedom. The 
asymptotic states thus feature a free electron with momentum p and energy e p — p 2 /2 (in atomic units) and an 
./V-electron molecular state \^q) \4>n) with the nth excited vibration or rotation in the electronic ground state with 
energy Eo n . For definiteness we choose the electronic kinetic energy as the zeroth-order electronic Hamiltonian _ff ,o 
in this section, describing the motion of free electrons. Note that the scalar product of Eq. (|52|) implies integration 
over all coordinates of the N + 1 electrons in addition to the integration over the nuclear coordinates R. 

The S'-matrix can be expressed by the particle part of the inelastic Green's function of Eq. ( JToj ) in the momentum 
representation as follows from the definition ( |l4| ) and Eqs. (|5^) and (|53"|): 

S(p'm «- pn) = Irrn^ e t ^' +E °^ Eoo)t G l ™ p n]+ (t,t') e -*(=-+^o»-«oo)*'. (54) 

t'—>— oo 

The particle part (j[ m <™]+ may be replaced here by the full Green's function Q^" 1 ' 1 ^ for two reasons: 

• The hole part Q^ m ^~ does not contribute because the S-matrix only contains information about the asymptotic 
region of scattering where Q^ m ' n ^ vanishes because a p \^^} vanishes. For the effective single-particle scattering 
equations we are going to derive in the following, this means that ensuring the correct outgoing boundary 
condition avoids contamination of the wavefunction by unphysical contributions. 

• The hole part Q^ n ' n ^~ does not contribute in the considered time limit because of the theta functions that appear 
explicitly in Eqs. (|l4|) and (|T6|). Therefore the particle part may be replaced by the full Green's function when 
the time ordering is treated properly. 
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We now introduce the inelastic improper self energy T(uS) in analogy to the common improper or reducible self 
energy [|l2]]|] by rewriting Eq. (|50|) as usual by 

£( w ) = £(°) (w) + £ (0) (w) 2» £ (0) (w), (55) 

2» = £{") + £ (0) M 2» ■ ( 56 ) 

Following the derivations of Refs. H|4l|l we can express the scattering S'-matrix by the improper self energy: 

S{p'm <- pn) = 8 nm 8 ppl - 2niT p [ 7 p n] {e p + E Qn ) 8{e p , + E Qm - e p - E 0n ). (57) 

The inelastic improper self energy 7~[ m, ™l at the scattering energy e p + E 0n can thus be identified with the on-shell 
T-matrix of scattering theory. Following standard scattering theory, we see from Eq. (^) that the inelastic self energy 
A^ m ' n \u>) (or equivalently the dynamical self energy A(u))) presents the (optical) potential for the scattering process. 
j\\m,n] j g -j-^g ptj ca i potential for multichannel scattering where the channels are defined by the eigenstates \<j> n ) of 
the nuclear motion. 

The relation of the inelastic improper self energy T to the improper self energy T of the usual (electronic) single- 
particle Green's function ( 0,0], see also Sj for explicit expressions that can be useful in the present context) can 
be described, in analogy to Eq. (EJ), by 

T} m ' n] (oj) = (ct> m \Z{u - H n + E m )\cp n ). (58) 

Note that this formal expression in fact expresses a rather complicated relation of the purely electronic improper self 
energy T(lo) and the inelastic scattering T- matrix because the nuclear kinetic energy in H n formally appears as an 
energy variable of T(w), which has itself a complicated dependence on the nuclear coordinates R. If, however, T n is 
omitted in this expression replacing H n by Vo(R), we can recover the well-known adiabatic-nuclei approximation (an) 
of electron- molecule scattering: 

T K„](a„) H = (0m | Tp/p(w _ Vq{R) + Em) | 0n) . (59) 



B. An effective Schrodinger equation 

In the preceding section we identified the inelastic self energy ^4l m ' n l with the optical potential for multi-channel 
scattering. Consequently, effective one-particle equations of the Lippmann-Schwinger or Schrodinger type can be 
derived. In the multichannel picture one obtains sets of coupled equations, one for each combination of initial and 
final channel. Instead of rederiving these equations, which can be found in Ref. [[ll| for the analogous case of 
electronically inelastic scattering, we will discuss directly the time-independent Schrodinger equation in the position- 
space representation. 

Any effective one-particle equation must involve explicitly the coordinates of the scattering electron and the labels 
of the nuclear states or, equivalently, the nuclear coordinates R. Instead of the full scattering wavefunction 

\*?{R)) = (R\% + )) 

= lim e-^+^-^at !<></>„(#), (60) 

t— > — oo 

which is the nuclear-coordinate representation of the scattering wavefunction of Eq. ( |53| ) and presents a wavefunction 
in the nuclear coordinates and an abstract ket in the N + 1-particle space, we introduce the effective or optical 
wavefunction f(r, R) by 

4" 1+ (r,H) = «|V(r) \*; + (R)). (61) 
This function can be calculated from the dynamical Green's function Q of Eq. ([l3]) by 

fW + (r,R)= lim [ dr'(R\G(r,0;r',t') |0 n ) <p p (r')e~^ P+ E 0n -E oo) ^ (62) 
where <p P (r) is a plane wave. The relation to the scattering S'-matrix is given by |5|,f4l[ 
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S(p'm ^pn) = JdR J dr <j>* m {R) ^,(r) /M+(r, H). (63) 

An effective Schrodinger equation for the considered scattering problem reads 

[H c .o + H n + A[E) - E 0Q ] /N+( Pj R) = E /M+(r, H), (64) 
as follows by analogy from Refs. [EllIM. A self-contained derivation of this equation using projection operators will 



also be given in Sec. |VC]. In the effective Schrodinger equation (|64j), E denotes the total energy with respect to 
the zero point defined at the molecular ground-state energy Eqq. This zero point is related to a situation where the 
scattering electron is so far away that it does not feel any forces exerted by the target, which is in its ground state, 
and the particle as well as the projectile are at rest. The zeroth-order electronic Hamiltonian H G $ is an operator 



acting only on the coordinates r of the scattering electron. According to the choice of Sec. IV A, it consists of the 
electronic kinetic energy and describes the motion of a free electron. Eq. (M), however, is also valid if the Hartree-Fock 
operator is chosen as the zeroth order. In this case, iJ e ,o contains also the static-exchange potential of the target 
molecule. The nuclear Hamiltonian H n acts only on the nuclear coordinates R and describes the nuclear motion 
of the target molecule in its electronic-ground-state potential Vo(R). The dynamical self energy A(E) of Eqs. (|i^ ) 
and ( J44] ) appears here in the coordinate-space representation and takes account of the complex many-body nature of 
the molecular target. It acts as a nonlocal (integral) operator in the electron coordinates r and also in the nuclear 
coordinates R: 

A(E) f(r, R) = J dR' J dr' A{r, R, r' , R' , E) f(r', R'). (65) 

Introducing the operator 

C(E) = H cfi + H n + A(E) - E 00 , (66) 

which can be seen as an analogue of the Layzer operator jl9) of the usual single-particle Green's function, the effective 
Schrodinger equation (Q) can be written as 

C(E)f(r,R)=Ef(r,R). (67) 

This is formally a pseudo-eigenvalue equation because the eigenvalue E also determines the operator C(E). In the 
scattering regime (E > 0) there is, of course, a solution for every E and the problem of solving Eq. (p7ft is to find 
a wavefunction f(r, R) with the correct boundary conditions such that the equation is fulfilled at a given scattering 
energy E. We usually identify the dynamical self energy A(E) with the dynamical optical potential although A(E) 
is not the only non-kinetic-energy component of C(E). 

As it stands, the effective Schrodinger equation ([54]) or (|67]) is an exact equation but, of course, approximations have 
to be introduced in most cases to determine the integro-differential operator C(E). We will discuss some possibilities 



in Sec. pV D[ The best strategy to solve Eq. (|64j) then certainly depends on the level of approximation introduced 
in the determination of C(E) and on the nature of the problem to be solved. A variety of methods is available. 
E. g. close-coupling p7| , p8[ or the other methods mentioned in Sec. | can be used to solve Eq. (^34|) directly for the 
scattering wavef unctions. Resonances may also be calculated using bound-state techniques b y e mploying the complex- 
rotation method ]5C|j5l|| , complex absorbing potentials [f^-jsil , or the stabilisation method [p5[ , etc. If necessary, the 
eigenvalues E may be found by iteration. Any other standard method for solving scattering problems will also do, 
since the corresponding Lippmann-Schwinger or effective time-dependent Schrodinger equation may also be used as 
starting points. 



C. A direct derivation of the scattering equations 



We will now show how to derive the effective Schrodinger equation (|6J) directly from the Schrodinger equation for 
the full scattering problem, which comprises the mot ion o f N + 1 electrons and the nuclear motion. This derivation 
is fully equivalent to the derivation sketched in Sec. IV A. Although the dynamical Green's function and its Dyson 



equation do not appear explicitly, the direct derivation, which uses projection operators in a composite Hilbert space 
with particle and hole states, is firmly based on Green's function theory and has to be seen as the result of recent 
developments in this area Q,^3). While traditional many-body Green's function theory supplies powerful methods 
to compute the dynamical optical potential approximatively (see Sec. IV D), the projection operator formulation of 
Dyson's equation j43|] yields a complementary way of reflecting on the Green's function approach and its relation to 
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Feshbach projection. Reconsidering in this section the algebra of Sec. Ill will help us to further clarify the role of an 
H-dependence in the electronic single-particle basis and in the electronic ground-state \^q)- We will also have an 
opportunity to discuss the role of unphysical contaminations of the optical wavefunction. 

We start from the time-independent Schrodinger equation for the full scattering problem of N + 1 electrons and 
the atomic nuclei: 

(H - E 00 ) |* to tGR)> - E \% ot (R)). (68) 

The total wavefunction l^totf-R)) is an arbitrary solution of the Schrodinger equation like, e. g., the scattering 
wavefunction of Eq. (|60|). The eigenvalue E measures the energy of the scattering system relative to the 

molecular ground state energy Eoq. Following Eq. (^), we split the total Hamiltonian H into the electronic part H e 
and the nuclear kinetic energy T n . Adding and subtracting the ground-state potential Vq(R) and introducing the 
nuclear Hamiltonian H n = T n + Vo (R) of Eq. (|Io|) , we may write 

[H e - V (R) + H n - E 0Q ] |*tot(fl)> = E |*tot(H)>. (69) 

By introducing the 2x2 matrix operator H of Eq. (^) we may rewrite Eq. ( p^ ) into the two-component equation 



where !L is the 2x2 unit matrix. Since the second component of the vector ( |^tot(-R)), 0)' is set to zero, Eq. (70) is 



completely equivalent to the initial Schrodinger equation (|68|). However, the vector ( \^ t ot(R)), 0) may also be seen 
as an element of the composite Y-space of Note |4q| . We now introduce the projection operator 

projecting onto the primary space, which is spanned by the vectors of Eq. ([23]). By definition, the projector P 

primarily acts on the electronic degrees of freedom. However, we also have to consider the dependence on the nuclear 
coordinates R because the nuclear kinetic energy T n contained in the nuclear Hamiltonian H n introduces derivatives 
with respect to R into Eq. ( |ro"| ) . The nuclear-coordinate dependence of P derives directly from the electronic ground 
state |*o), 

as can be seen from the definition (^3|) of \Y p ) when we assume that the basis of single-particle orbitals 
tp p (r) is independent of R. Since the primary space and consequently the projector P are invariant under a change of 
the single-particle basis, this also holds true when U-dependent orbitals like, e. g., Hartree-Fock orbitals are chosen, as 
long as they from a complete basis of the single-particle space. Consequently a formulation employing the projection 
operators P and Q = 1 — P remains independent of the choice of the single-particle basis and thus an independence 
of the orbitals does not matter. 



We already discussed in Sec. HI B that we want to adopt the Born-Oppenheimer approximation for the electronic 
ground state l*^}. In particular we assume the Independence of the electronic ground state to be weak. As 

a consequence, the projection operator P commutes with the nuclear kinetic energy T n and therefore also with H n : 

[P,T n ] = [P,H a }=0. (72) 

We want to stress once more that we assume the Born-Oppenheimer approximation only for the nuclear motion in 
the electronic ground state of the A^-electron molecule, which is physically reasonable for many molecules. However, 
this neither implies adiabaticity for the scattering electron nor for excited states of the scattering complex. By way 
of contrast, the nonadiabatic coupling of the projectile motion to the nuclear motion, which can be very important 
for slow projectiles, is explicitly accounted for in our formalism. 

Defining the operator Q = 1 — P as a projector onto the secondary space (1 is here the identity operator in Y-space 
p5j), we can easily partition Eq. (|7^) in order to obtain an effective equation in the primary space. We take the 
following steps: 

Inserting the identity 1 = P + Q into Eq. (|7^) yields 

HP + Q) ( l * tOt m )=(E-H n + E 00 ) ( ) . (73) 

Acting on this equation with Q from the left and using [Q, H n ] = 0, which follows from [P, H n ] = 0, yields an equation 
for the secondary-space component of the total wavefunction Q ( \ ^ t ot(R)), 0) : 



14 



HqpP 



|*tot(fl)> 




(E - H n + E 00 - Hqq) Q 



|*tot(fl)> 





(74) 



Here we introduced the notation Hqp for the operator product QHP, etc. Using Eq. (|7J) to replace the secondary- 
space component from the P-projection of Eq. (73) leads to the desired equation for the primary-space component of 
the wavefunction: 



Hpp + H n — Eqo — Hpq 



1 



E-H„ 



En 



H 



-H, 



QP 



QQ 



l*tot(«)> 





= EP 



l*tot(fl)> 




(75) 



This is the resulting form of the time-independent Schrodinger equation projected to the primary space. The projected 
component fulfils a pseudo-eigenvalue equation with an energy-dependent operator on the left-hand side. This operator 
acts on the nuclear coordinates R and on the electronic coord inate s, but only within the primary space, which is 
isomorphous to the one-particle space (as was discussed in Sec. [II A ). 

In the following we will make contact to our previous formulation and identify the self energy of the dynamical 
Green's function as the optical potential. The primary-space component of the wavefunction is given by 



P 



l*tot(fl)> 





Y„ 



l*tot(H)) 




\Y q )f g (R), 



(76) 



where f q (R) = (^\a 
inner product of Eq, 



'i'tot(-R)) is the effective wavefunction of Eq. 
75) with (Y p \ from the left yields 



( |6l|) in the orbital representation. Taking the 



({Y P \H\Y q ) + (H n -E 00 )5 p 



(Y p \HQ 



E — H n + Et 



00 



He 



-QH\Y q ) f q (R) = Ef p (R) 



(77) 



The first term on the left can be identified as the primary block H_ qq of Eq. (]2§|) , which splits into the zeroth order 
energy matrix e and the static (i. e. energy-independent) part of the self energy -4(oo) according to Eq. (^3|). Also 
the energy-dependent part of A(u>) can be identified in Eq. ( \n\) when the representation 



(78) 



is chosen for the projector to the secondary space. The sum runs over the basis states \Qj) of the secondary space 
only [c. f. Eqs. @) and ©]. We obtain 



J2 i £ Pi + ( H » ~ £ oo)5 M + A pq {uo)\ f q (R) = Efp(R). 



(79) 



The effective Schrodinger equation in the form of Eq. (|64|) is obtained when using the coordinate representation instead 
of the orbital representation for the scattering electron and realising that the matrix e is the orbital representation 
of the operator H e n of Eq. ( |32| ) in the one-electron space. As discussed earlier in this subsection, the projection 
operators P and Q are completely independent of the choice of the single-particle basis and thus Eq. (^9|) holds in the 
momentum, in the coordinate-space, and in arbitrary orbital representations. We are also not restricted to a particular 
choice of the zeroth-order electronic Hamiltonian i? o ,0 ; represented by the matrix s in Eq. (|79|). Any single-particle 
operator, which may or may not depend on R, is possible. A convenient choice for the electron-molecule scattering 
problem is certainly to choose the Hartree-Fock operator for H e ^ because this simplifies the calculation of the self 
energy A(ui). In this case, H c n, contains the electronic kinetic energy as well as the static exchange potential. 

So far, we have proven that the projection f(r, R) of Eq. (|6l]) of the physical scattering wavefunction fulfils the 
effective Schrodinger equation (p4j) . There are also unphysical solutions of this equation due to the introduction of 
the second component in Eq. fl70|). The unphysical component lives in the Hilbert space of N — 1 electrons, coupled 
to the same nuclear degrees of freedom as the physical component. Since the Hamiltonian for the second component 
from Eq. (70) is given by 



V (R)-H e + H n -E l 



oo 



-[H e -T n -2V (R)+E 00 ], 



(80) 



the nuclear dynamics is not treated in the correct way to describe an ionised molecule because H e and T n have 
opposite sign. If we want to apply the formalism of the dynamical optical potential to describe the dynamics of 
ionised molecules, the definition of the Hamiltonian has to be changed accordingly. 
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Coming back to the description of the electron-scattering problem we want to consider the question of how to 
identify unphysical solutions and whether the unphysical component can interfere with the physical solutions. The 
answer to these questions is easily given when the scattering boundary conditions are obeyed correctly: Since the 
electronic ground state of the scattering target l^^} represents a bound state, any overlap (^| ip(r) I^q) is a square- 
integrable function of r vanishing asymptotically for |r| — > oo. Consequently, the same holds true for the effective 
wavefunction / un P h y s (r, R) belonging to an unphysical solution of Eq. (70) 



/ 



unphys 



(r,R) 



Y(r) 





N-l 
unphys 



OR) I 



N-l 



(81) 



as a function of r. As long as one solves the effective Schrodinger equation (|64|) or the corresponding Lippmann- 
Schwinger equation imposing scattering boundary conditions on the effective wavefunction f(r,R), one obtains a 
physical solution. Even a contamination of the physical solution with a short-range unphysical solution will do no 
harm because for the calculation of the S'-matrix elements only the asymptotic behaviour of the effective wavefunction 
enters. 

Let us briefly examine the only approximation made in the derivation of Eq. (|T5|). Without the approximation (|7^), 
which allowed the commutation of the nuclear kinetic energy T n with the projection operators P and Q, additional 
terms proportional to PH n Q and QH n P appear in the projected Schrodinger equation. These terms describe virtual 
excitations from the primary to the secondary space, which mean an electronic excitation in the target molecule, 
mediated by the nuclear kinetic energy operator. These terms compete with excitations by the electronic Hamiltonian 
PH C Q and QH C P against which they can usually be neglected. The only exemption may be given in cases where 
the Born-Oppenheimer approximation for the electronic ground state fails and the potential surfaces of the electronic 
ground and excited states of the target molecule come close to each other. In these rare cases, the vibronic interactions 
have to be accounted for and a diabatic representation for the ground state may be adequate. The dynamical optical 
potential then has to be augmented by terms that describe the corresponding vibronic transitions. 



D. Approximative scattering equations 

In this subsection we want to discuss different approximations to the exact scattering equations and the physical 
models these approximations imply. In the first two points we will consider possibilities to regain well-known approx- 
imate descriptions of electron-molecule scattering by either neglecting the energy-dependent optical potential or by 
freezing the nuclear degrees of freedom in the exact equation (p7|). In the last two points we will concentrate on the 
dynamical optical potential itself. Apart from discussing the role of the nuclear dynamics within the optical potential 
we will show how to compute ab initio dynamical optical potentials. 



1. Static approximations to the optical potential 

First of all we want to consider the simplest of all approximations to the optical potential, which results by neglecting 
the optical potential A(lo) altogether. Formally this approximation is equivalent to neglecting the interaction part 
fle.i in the electronic Hamiltonian H e of Eq. (fUj). This equivalence follows easily from Eq. (^4|) which connects the 
dynamical self energy A(uj) to the usual self energy which is by definition at least of first order in H c , \ [n2l . 

In other words we can say that, neglecting the dynamical self energy A(lu), the effective Schrodinger equation (|6j)is 
still exact for physical systems where the scattering electron is not correlated with other electrons in the system and 
thus is well described by iJ c .o- The one-particle operator H c q may still contain mean fields or forces exerted by the 
atomic nuclei or external fields, of course. 

The level of approximation gained with the resulting static (i. e. energy-independent) Layzer operator C st = 
H c q + H n — Eqq depends on the choice of the zeroth-order electronic Hamiltonian H c : 

• Choosing H c ,o to describe the electronic kinetic energy yields a Layzer operator C st that describes the separable 
motion of a free electron and a vibrating or rotating target molecule in its electronic ground state. 

• Choosing instead H c $ to describe the independent Hartree-Fock operator yields a static Layzer operator that 
describes the motion of the scattering electron under the influence of the static and the exchange interactions 
with the Hartree-Fock charge cloud of the target molecule coupled to the nuclear motion in the usual electronic- 
ground-state potential Vo(R) augmented by the Coulomb attraction between the nuclei and the electronic 
projectile. This approximation gives a consistent treatment of the electron-molecule scattering problem in the 
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strict single-particle picture. Projectile and nuclear motion are fully coupled but the very nature of this static 
approximation excludes polarisation effects as well as the possibility of electronic excitation of the target caused 
by the impact of the scattering projectile. Another defect of this Hartree-Fock based static-exchange approxi- 
mation is that the static charge cloud of the target molecule is only described in an approximate (uncorrelated) 
manner through the Hartree-Fock wavefunction. 

Another straightforward possibility for an energy-independent approximation to the dynamical optical potential 
is to include the high-energy limit of the self energy A(oo), which is also called the static self energy. The static 
part of the dynamical self energy A(oo) is identical to the purely-electronic static self energy E(oo), as can be seen 
from Eqs. ( |35| ) and ( [f3| ) [see also Eq. (fli|)]. As shown in Ref. [^6), this nonlocal operator can be interpreted to 
improve the static-exchange interaction beyond the Hartree-Fock description of the target wavefunction including 
target correlation. The present approximation to the Layzer operator C csc = H c q + H n + A(oo) — Eqq is independent 
of the choice of the zeroth-order electronic Hamiltonian H e jj in Eq. ( |3l| ) because the static self energy takes account 
of the interaction term ff e ,i- The Layzer operator £ csc corresponds to a scattering potential that originates from a 
correlated but static charge distribution of the target molecule, which is also known as the correlated static exchange 
(cse) potential. 

Approximations to the energy-dependent pa rt of th e dyn amical optical potential in the context of the fully coupled 



scattering problem will be discussed in Sees. IV D 3 and IV D 4. In the following paragraph we will consider the 



simplified case where the nuclear dynamics is neglected altogether. 



2. Fixed-nuclei scattering 

Let us recall that the full Layzer operator C may be expressed with the help of Eq. ([lij) by 

C(E) = H e , + T n + V Q (R) + £(£ - T n - V (R) + E Q0 ) - E m . (82) 

In the limit of infinitely heavy nuclei, the nuclear kinetic energy T n can be neglected. We call this limit the fixed-nuclei 
(fn) limit because the Layzer operator 

C hl (E) = H c , + E(E- V (R) + E 00 ) + V (R) - E 00 (83) 

and the corresponding scattering wavefunction /(r,i?) of Eq. (Q) now depend only parametrically on the nuclear 
coordinates R. The operator C lu (E) is exactly the Layzer operator used in the literature for electron scattering from 
atoms or rigid molecules fll|,ff9|]5|] apart from the Independent energy shift Vq(R) — E o- This energy shift resets the 
zero point of the energy scale from Ego to Vo(R). 

In the present case of fixed-nuclei scattering, the optical potential is given by the purely electronic self energy S(w). 
As mentioned earlier, the static part £(oo) improves the static-exchange potential with respect to the correlation of 
the t arg ets electronic ground state wavefunction. The energy-dependent part M(u>) — — E(oo) is, according to 
Eq- @, g iven b Y 

M(uj) = H u l -^H u . (84) 

— =a ojI- H_ bb 

For convenience, we again resort to the matrix notation for electronic coordinates used already in the preceding 
chapters. Of course, the matrices can also be expressed in the coordinate representation. M(uj) accounts for the 



so-called dynamic correlation including the polarisation of the target by the incident projectile electron 57 1. The 



energy dependence of M(u) has been studied in the literature J 58| and the consequences of the energy dependence of 



optical potentials for scattering systems have been discussed 



Here, we briefly want to remark the following: 



M{uj) introduces poles and branch cuts in the optical potential, which can be found by diagonalising the denominator 
of Eq. (p4]). The matrix H_ bb in the denominator is already diagonal in zeroth order, e. g. in the Hartree-Fock 
approximation. The diagonal elements correspond to electronic excitations of the (N + 1)- and (N — l)-electron 
systems and can be classified by configurations of 2p — h, 3p — 2h, etc. character as well as configurations of p — 2/i, 
2p — 3h type, and so on. Explicit matrix representations for H , H , and H_ bb can be found in Ref. |42]| . 

For the calculation of S(w), standard approximation schemes are available. The so-called algebraic diagrammatic 
construction (ADC) scheme [ H , for example, yields a perturbation-theory-based hierarchy of approximations for 
E(w) that preserve the analytic structure of this function of u>, in contrast to ordinary perturbation theory, which 
does not. The ADC approximation to S(w) has the structure of Eq. (|35| ) but the matrices H are constructed in 
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the nth order ADC scheme such that the perturbation expansion of the approximate E ADC (a->) with respect to the 
electronic interaction iJ .i of Eq. (|3l| ) coincides up to nth order with the usual Feynman-Dyson perturbation series 
ofE(w). 

Since the poles and cuts relating to the ionised molecule (N — 1 electrons) appear at negative energies, they only 
have a weak influence on the energy-dependence of the optical potential at the relevant scattering energies. The poles 
and cuts at positive energies, however, can be associated with excited states of the (N + l)-electron scattering system. 
Signatures of Feshbach resonances therefore can be found in the dynamic part of the optical potential although the 
exact resonance position and width, of course, have to be calculated from the full effective Hamiltonian. A branch 
cut will appear in the analytic structure of above the first excitation energy of the target molecule because with 
sufficient energy for an electronic excitation, a new channel opens and inelastic scattering becomes possible. The 
branch cut of E(u>) has the consequence that the optical potential acquires an imaginary component and becomes 
non-hermitian, thereby accounting for the loss of amplitude in the channel of electronically elastic scattering. 



3. Nuclear dynamics in the optical potential: Expansion with respect to T n 



We are now going to discuss the influences of the nuclear motion on the dynamical optical potential A(E). The 
dynamical self energy A(E) as given by Eq. ( fi3| ) consists of an energy-independent (static) part A(oo) and an energy- 
dependent part A4(E) = A(E) — _4(oo), which vanishes for large energies E. The energy-dependent part A4(E) 
includes all effects of dynamic correlation including polarisation but is now modified with respect to the fixed-nuclei 

the energy-dependent part reads in matrix notation: 



case considered in Sec. IVD 2| . As follows from Eq. 



M{E) = H 



1 



[E-T n -V (R) + E 00 ]l-R 



(85) 



Like before, T n denotes the nuclear kinetic energy, Vq(R) the (electronic) ground-state potential, E is the scattering 
energy, and E 00 the molecular ground-state energy. The matrices -ff . are the partial matrix representations of the 



two-component electronic excitation-energy operator H, like introduced in Sec. Ill A , and depend on the nuclear 
coordinates R. In order to diagonalise the denominator to find the pole structure of Ai(E), not only eigenvectors of 
H hh corresponding to electronic eigenstates in the secondary space are necessary but instead one also has to consider 
the additional degrees of freedom of nuclear motion. Since the energies usually associated to the nuclear dynamics are 
much smaller than typical energies of electronic excitations, we may assume that th e coarse structure of the energy 
dependence does not differ much from the fixed- nuclei case discussed in Sec. IV D 2 However, the nuclear dynamics 
will introduce a fine structure and may become very important, either at energies where electronic excitation is 
possible or nearly possible, or in situation that are very sensitive to disturbance of the scattering potential, e. g. close 
to vibrationally inelastic thresholds. 

As long as the scattering energies are safely away from Feshbach resonances or inelastic thresholds, the influences of 
nuclear dynamics on the optical potential itself may either be neglected or accounted for approximately by expansions 
of the denominator in Eq. (|85|) . We will now proceed to discuss two convenient expansions for M(E) . 

Closest in spirit to the before-discussed fixed-nuclei approximation is the assumption that T n represents a small 
perturbation, suggesting the following type of expansion: 



M{E) = H n 



1 



E K 



i 



= ab [E - Vo(R) + E Q0 }1 - R bb ^ \"[E- V (R) + E 00 }1 - K bb , 
In the present context, we call the first term of this expansion (y = 0) the adiabatic optical potential (aop): 



(86) 



A aop (E) = H -e + H , 



i 



[E-Vo{R) + E 00 }l-K bb 
UE-V (R)+E 00 ). 



(87) 



For mally, th e adiabatic optical potential A aop (E) coincides with the optical potential used in fixed- nuclei scattering 



in Sec. [VP 2. The difference is that the nuclear coordinates R are now dynamical variables. The adiabatic optical 



potential „4 aop (i?) is a local operator on the nuclear coordinates R because Vq and the matrices H depend on R 
but there are no derivatives with respect to R in A aop (E). Since A aop (E) can be calculated, for each R, from the 
purely electronic s elf ener gy E(i£), standard techniques for calculating the electronic self energy like the ADC method 
mentioned in Sec. [VP 2 can be applied. 
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The Layzer operator corresponding to the presently discussed approximative treatment of the dynamical optical 
potential „4 aop (.E) is given by 

£ aop (£) = H cfi + H n + A aop (E) - E 00 . (88) 

Although the dynamical optical potential is treated here in an adiabatic approximation, the Layzer operator C aop (E) 
describes the fully and non-adiabatically coupled motion of the projectile electron and the atomic nuclei with a 
simplified optical potential. In this approximation, the nuclear dynamics in the electronically excited states of the 
target is treated adiabatically. Even though this approximation may not be justified when describing experiments 
probing in detail the electronic excitation structure of the target, the situation is different for scattering at low but not 
too low energies. When electronic excitations are forbidden by energy conservation because the electron-impact energy 
is too low, the adiabatic approximation „4 aop (_E) may be very good. At low impact energies, electronic excitation then 
is still possible as a virtual process leading to polarisation of the molecule but obeying the energy-time uncertainty 
relation. When the typical times for these virtual excitations are much smaller than the time scales of nuclear 
vibrations, which will usually be the case, then the fixed-nuclei approximation for the dynamical optical potential 
„4 aop (£;) is adequate. 

For very slow scattering electrons, on the other hand, it may also be necessary to improve the approximation 
A aop (E), when the cross sections probe the details of the optical potential very sensitively. In the adiabatic optical 
potential A aap (E), the dynamical relaxation of the nuclear structure during the polarisation of the electronic charge 
cloud induced by the incoming projectile is not treated properly. At very low projectile velocities where the trajectories 
are most sensitive on the coupling to the nuclear motion this may become an important defect of the theory. In this 
case the expansion of the dynamical optical potential discussed in the next subsection, which allows for an approximate 
treatment of the neglected effects, becomes most valuable. 



4- Expansion of the dynamical optical potential with respect to H n — Eqo 



The nuclear kinetic energy operator T n introduces derivatives with respect to the nuclear coordinates R in the terms 
of the expansion ( [B6|) for v > 1. In general, these terms will be difficult to evaluate because the approximate and 

also the exact matrices H_ . . may depend strongly on the nuclear coordinates R due to the electron- nucleon Coulomb 
— i j 

repulsion contained in H. Another expansion that is better suited for higher-order approximations is obtained when 
taking T n + Vo(R) — Eqq = H d —Eqo as the "small" perturbation. The assumption that H u — Eqq may be regarded small 
in the denominator of Eq. (|85| ) is justified when the energies associated to nuclear excitations appearing during the 
scattering process and possibly during (virtual) electronic excitations are small compared to the electronic excitation 
energies of the target. The corresponding expansion of the energy-dependent part of the dynamical self energy reads 



M(E) = H 



1 



E 



nb El-H^ 

= =bb f=0 



(Hn — E 



1 



H, 



=ba 



(89) 



The first term Ai (E) of this expansion [y — 0) can be identified as the energy-dependent (dynamic) part of the 
purely electronic self energy ( |8~i| ) and will be called the zero-point optical potential because it reflects the optical 
potential of ground-state nuclear motion: 



M°(E) = H 



1 



K ba =M(E) 



ab El-H_ bb 

The zero-point optical potential is a local operator with respect to the nuclear coordinates R like the adiabatic optical 
potential A4 aop (E), which was discussed above. 

The higher order terms in the expansion for v > 1, however, are differential operators in the nuclear coordinates. 
They can be transformed into nonlocal integral operators M"(E) f(R) = J dR' M U (E, R, R') f(R') by inserting the 
resolution of the identity with respect to nuclear motion Y^k~~\ ( t ) k) (4>k \ ■ In the coordinate representation for the nuclear 
degrees of freedom R, the integral kernel belonging to the second term of the expansion (pj) reads 



M}(e,r,r') = J2kJR) 



E l-iLJ R 



: (R)[E { 



Ok 



E i 



on J 



KM 



l\=ba 



(90) 



where we explicitely indicated the nuclear-coordinate dependence of the matrices H_ . . The numerical realisation of this 
approximation is not very difficult because the nuclear eigenvalues i?ofc and functions tfik (R) can be calculated easily 
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form the ground-state potential Vq(R) and approximate matrices H are available, e. g., by the ADC approximation 



as described above in Sec. 1VD 2 



Eq. (J90j) is remarkable because it allows us to estimate the importance of the higher-order terms in the expansion 
(|S9|). If the nuclear dynamics of the target molecule, in spite of virtual excitations by the scattering electron, can 
be approximated well by the zero-point vibration of the isolated molecule (in its electronic ground state) then all 
terms with k ^ in Eq. ( |90| ) can be neglected. Since, however, the energy difference [E k — E 00 ] vanishes for k = 0, 
the term M}{E) and the higher order terms in ( |S9| ) vanish. In this approximation only the zeroth term M°(E) 
contributes in the expansion of the dynamical optical potential. In other words M_{E) is a local operator in the 
nuclear coordinates if the nuclear-coordinate dependence of the effective scattering wavefunction f(r,R) is given by 
the nuclear eigenfunction (f>o(R) and the R dependence of the matrices £i bb and H_ b can be neglected. This will not 
always be the case. However, it can also be tested during a calculation for which values of k the integral 

1 

L 



dRMK) — -g — R ba (R)l(R) (91) 



gives the largest contribution. The approximation M. (E) can then be improved by expanding around H n — E, where 
E is a suitably chosen mean vibrational (rotational) energy, instead of H n — Eqq, yielding the zeroth approximation 

M° = £(£ + E q — E) — £(co) (92) 

and the corresponding higher approximations like in Eq. (|8S|). 

The difference between the adiabatic optical potential A aop (E) and the zero-point optical potential A°(E) or 
the improved approximation A (E) lies in the energy reference of the denominator of the optical potential. We 
expect A°(E) or A°(E) to be suitable approximations for A{E) when the nuclear motion is more or less confined 
to the immediate surroundings of the equilibrium configuration H eq of the molecule. In particular for processes 
that travel between different levels on the potential energy surface Vo(-R) corresponding to a considerable transfer 
of energy between the projectile and target like it happens in associative detachment or dissociative attachment, 
we expect the adiabatic approximation A aop (E) to be the better-suited approximation. However, we want to stress 
once again that the treatment of the nuclear dynamics within the optical potential is certainly less essential than 
the direct coupling between projectile and nuclear motion in the effective Schrodinger equation (p4|), as long as the 
scattering energy remains well below the threshold for electronic excitation of the molecular target. The expansion 
( p9| ) presents a systematic possibility to improve upon the adiabatic or the zero-point approximation that is applicable 
if the scattering process has too little energy to access the electronically excited states of the target. The first-order 



term A4 1 of Eq. (90) and the corresponding higher approximations provide access to a systematic improvement of 
the discussed approximations for ab initio methods of quantum chemistry that may prove necessary to describe the 
dynamical couplings correctly for very-low-energy scattering processes that are particularly sensitive to the correct 
description of the optical potential. If the scattering energy is above or close to the threshold for electronic excitation 
of the target molecule, the effective Schrodinger equation ( p4| ) still remains valid but other approximations to the 
dynamical optical potential have to be used. 



V. CONCLUSIONS 



In this paper we have derived a rigorous optical potential for the coupled motion of the projectile electron and 
the atomic nuclei in electron-molecule scattering. The dynamical optical potential should be particularly valuable for 
studying the non-adiabatic coupling of projectile and nuclear motion close to inelastic thresholds where the projectile 
velocity is comparable to the typical velocities of nuclear motion. Our work extends the well-known many-body 
formalism of optical potentials beyond the fixed-nuclei approximation. On the other hand the rigorous derivation 
and the explicit expressions given for the dynamical optical potential provide a justification and prospective for 
enhancement of model polarisation potentials used commonly in close-coupling and i?-matrix calculations. 

The dynamical optical potential can describe all kinds of electronically elastic scattering processes such as elastic 
scattering, vibrational or rotational excitation (de-excitation), dissociative attachment, and associative detachment. 
Also when electronic excitations are energetically accessible, the dynamical optical potential may be used to calculate 
the (electronically) elastic partial cross sections as well as resonance positions and widths. 

One central result of the present work is that for scattering energies far enough away from electronic excitations of 
the target, the dynamical optical potential is given in a first approximation by the static- exchange potential augmented 
by the fixed-nuclei self energy £(w) of the purely electronic Green's function. This approximation, which yields a 
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local operator with respect to the nuclear coordinates R, can be improved by a hierarchy of non-local terms. Due to 
the close relation to the usual self energy all terms can be calculated with standard ab-initio methods using 

well-defined approximation schemes. 

The dynamical optical potential provides a firm basis for theoretical studies of coupled projectile and nuclear motion 
in ultra-low-energy electron-molecule scattering, which is becoming a vivid field of interest due to recent experimental 
advances (see e. g. JH^]). The effective scattering equations presented in this paper can be solved, e. g. by using 
close-coupling expansions or other time-independent techniques. We expect that particularly valuable insight into the 
mechanisms of the scattering process, however, can be gained by time-dependent wave-packet calculations. Solving 
these scattering equations certainly remains a numerically demanding task but becomes more and more feasible now 
with state of the art computing facilities and the advancement of numerical methods |6C| l . 
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